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Abstract 

We study cosmological evolution after inflation in models with non-minimal deriva¬ 
tive coupling to gravity. The background dynamics is solved and particle production 
associated with rapidly oscillating Hubble parameter is studied in detail. In addition, 
production of gravitons through the non-minimal derivative coupling with the inflaton 
is studied. We also find that the sound speed squared of the scalar perturbation oscil¬ 
lates between positive and negative values when the non-minimal derivative coupling 
dominates over the minimal kinetic term. This may lead to an instability of this model. 
We point out that the particle production rates are the same as those in the Einstein 
gravity with the minimal kinetic term, if we require the sound speed squared is positive 
definite. 
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1 Introduction 


Inflation theories, now accepted as a successful way of explaining the isotropy and homo¬ 
geneity as well as the density fluctuation of the universe, assume the existence of a scalar 
held, or inhaton, to trigger inhation. The coupling between the inhaton and gravity might 
in general be non-canonical, and among such theories there exists a class which introduces 
no extra degrees of freedom. This class of models has been investigated especially after the 
discovery of Higgs boson at LHC nua and after the results from Planck satellite [3], in the 
context of having successful inhation with the standard model (SM) Higgs held. For exam¬ 
ple, it was proposed to couple the Higgs held non-minimally to the Ricci scalar R iHS] , or to 
introduce a non-minimal derivative coupling between the Higgs held and the Einstein tensor 
G,. IT]. Moreover, it was pointed out that a modihcation of the Higgs kinetic term from the 
canonical one also leads to a successful inhation [8]. For general setups for inhation using 
non-minimal couplings without introducing additional degrees of freedom, see the studies of 
the generalized G-inhation [9]. 

While inhation itself is a matter of great interest, the reheating phase is equally interesting 
not only because it necessarily follows the inhationary phase in order to produce SM particles, 
but also because it is expected to have rich phenomenology. To study the reheating phase, 
we must know the background evolution in the inhaton oscillation dominated universe and 
also how the inhaton couples to SM helds. In models with extended gravity, these aspects 
are often non-trivial compared with the standard inhation models with Einstein gravity. 
Interestingly, it is sometimes possible to reheat the universe without introducing explicit 
interaction terms among the inhaton and SM helds, since the inhaton necessarily couples to 
SM sector gravitationally at the tree level HQHn]. 

In this paper, we study the reheating in the model with non-minimal derivative coupling: 
C ~ G^^d^cjydycj). We do not necessarily assume that 0 is the SM Higgs held. The detailed 
analysis of inhaton oscillation dominated universe in this class of model was performed 
in Ref. [13] • There it was found that the Hubble parameter rapidly oscillates if the non- 
minimal kinetic term of the inhaton takes a dominant role. Thus it is expected that particle 
production of non-Weyl invariant light helds may be efficient. More interestingly, the inhaton 
has a large non-minimal coupling to graviton, which opens up a new possibility to produce 
signihcant amount of gravitons from the inhaton coherent oscillation. Actually we will hnd 
that gravitons may be resonantly produced after inhation. However, we will also point out 
that the sound speed squared of the scalar perturbation rapidly oscillates between positive 
and negative values in such a case and it may lead to an instability for the shortest wavelength 
mode of the scalar perturbation. 

The organization of this paper is as follows. First we write down background equations 
of motion, and summarize diherent phases the universe undergoes in this model. In Sec. E] 
we study the perturbation around the background and hnd that the sound speed squared 
of the scalar perturbation can be negative. In Sec. IH we analyze particle production in this 
model, focusing on minimally-coupled scalars and gravitons. Sec. |5]is devoted to conclusion. 
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2 Background dynamics 


2.1 Background equations 

First of all, we derive the background equation [7]. We consider the following action 


S= d'^XyZ—g 


2 2 




M2 


d^(j)d^(j) - V ( 0 ) 


( 2 . 1 ) 


where 0 is a real scalar held, or the inhaton, R is the Ricci scalar and G/j,^ = — g^^R/2 is 

the Einstein tensor. Also, Mp = is the reduced Planck mass and M is some mass 

parameter. We assume the background part of the metric to be the Friedmann-Lemaitre- 
Robertson-Walker (FLRW) one with negligible curvature 

ds^ = —dt^ + a{tYdx^dx\ ( 2 . 2 ) 


where a is the scale factor. After integration by parts, the background action becomes 


S 


d'^x oZ 





2 


R(0) 


(2.3) 


where M = a/a is the Hubble parameter. 

The Friedmann equation is obtained by introducing the lapse function N hy dt ^ Ndt 
and taking variation with respect to it: 


= 

P<f) = 


P<t> 

3Mj,' 




(2.4) 

(2.5) 


where is the energy density of cj) and we have taken A^ = 1 after the variation, 
the Friedmann equation 02.41) implies the following relation: 

_ 0^ 2 

^ M2M2 3‘ 

The variation with respect to the scale factor gives 


Note that 

( 2 . 6 ) 


3H‘^ + 2H 


M|.’ 






1 d 
JpJt 




(2.7) 

( 2 . 8 ) 
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where is the pressure of <p. As we will see later, violent oscillating features of 0’s energy- 
density and of the Hubble parameter for if S> M originate from the last term in Eq. fl2.8p . 
The equation of motion for 0 is 


1 + 


3H^ 


, 2H 


dV 


(2.9) 


From this equation of motion, we can see that the potential is effectively suppressed due to 
the non-minimal derivative coupling for H/M ^ 1. Note that Eq. fl2.9p is equivalent to the 
energy conservation equation 


P<j> + 3ii + p^) — 0, 


( 2 . 10 ) 


and that one of the three equations 02.41) . 02.7p and 02.91) is redundant. Also, H is calculated 
by eliminating 0 from the time derivative of Eq. 02.4p and Eq. 02. 9 p as 


H 


(1 + ^)(1 + 


HV 

>2^ 


(i + S-a 


9H^\ e ’ 
M2 ) 2 


where the prime denotes the derivative with respect to 0. 
In the following, we consider the power-law potential 

1 / = - 0 ”. 
n 


Here we define the effective mass rries as 



= f 


(f » 1) 
(f « 1) 


( 2 . 11 ) 


( 2 . 12 ) 


(2.13) 


where <I> is the amplitude of the 0 oscillation. Note that the effective mass is an increasing 
function with time due to the suppression factor M/H for H/M 1. Also note that rries^ 
remains constant for H/M ^ 1 in the inflaton oscillation dominated era. We sometimes 
express as for n = 2. 


2.2 Evolution of the universe 

The evolution of the universe for models with the non-minimal derivative coupling was 
investigated in detail in Ref. [I3]. Here, we briefly summarize the discussion given there. We 
assume that the universe is dominated by the inflaton 0 and neglect the damping of 0 due 
to particle production in this subsection. 

First, consider the case where meg <C H (Phase 0). In this case, inflation takes place due 
to the large Hubble friction term. We summarize the inflationary predictions of the models 
with non-minimal derivative coupling in App. |Al 
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Next, consider the case where M H rrieg (Phase 1). In this case, 0 starts to oscillate 
while the kinetic term is dominated by the non-minimal one. In analyzing this phase, we 
note that the energy density of 0 oscillates violently with time. To see this, one may refer 
to Eq. fl2.10p . which is given after substituting the explicit forms of p,j, and as 


A 

dt 




+ 3H 



3H d 

JpJt 



0 . 


(2.14) 


The second term is 0{Hp^), and no oscillating feature is caused by this term. On the other 
hand, the third term is 0{mef[P<f>), and hence p^ oscillates with the frequency of ~ meg due 
to this term. Note that it means that the Hubble parameter also violently oscillates with 
time when 0 dominates the universe. Because of this oscillating feature, does not suit for 
the evaluation of the evolution of the universe in this regime. Instead, it is shown in Ref. [13] 
that the following quantity J is an adiabatic invariant]^ 





3MlH 


1 \ 

2 M|,M2 j ■ 


(2.16) 


We prove the adiabaticity of J in App.jBl We can see from this equation that the oscillation of 
the Hubble parameter and that of 02 correlate with each other. In Ref. na, J is extensively 
used to calculate the expansion law of the universe in this regime. The expansion law in this 
regime was derived as 


{H) = 


2?7, + 2 1 


(2,16) 


3n t ’ 

where the parenthesis denotes the time average. We further recall the important feature in 
Phase 1 here. Notice that the Friedmann equation 02.41) implies that 


m 


eff 


$2 


Mj,M^ 


0 ( 1 ), 


(2,17) 


is always kept in Phase 1 when 0 dominates the universe. This is because the amplitude 
<h decreases while the effective mass rries increases with time, and hence rries^ is almost 
constant during Phase 1. This feature will become important when we discuss particle 
production in Sec. ID 

Finally, consider the case where H M (Phase 2). In this case, the kinetic term of 0 is 
dominated by the minimal one, and the evolution of the universe approaches to the ordinary 
case of the minimal oscillating scalar with Einstein gravity. The expansion law in this case 


is given bjl^ 


H = 


n 


3n 


21 

t 


(2.18) 


Throughout this paper we call a quantity I an adiabatic invariant if I satisfies I ~ 0{HI). 

In fact, for M{MM, the term proportional to 0 is still dominated by the term 
~ H /M^ in Eq. (12.91) . However, this term does not act as friction because it is not positive definite, and the 
expansion law soon becomes the ordinary one described by Eq. (12.1811 after a few Hubble time after entering 
the phase 2. 
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Whether Phase 1 exists or not depends on the model parameters. In order to discuss 
each situation separately, we consider the following two cases: 


Case A : Phase 1 does not exist. This corresponds to the case where H < M and 
meff < M are always satished after inflation. It requires 


AM"-2 

M2 


< 1 , 


(2.19) 


(see Eq. flA.9j) h Then, the cosmic evolution is the same as that of the minimal scalar 
held with Einstein gravity. Note that we always have H ~ in this case. 


• Case B : Phase 1 exists. In this case, the inequality fl2.19p is inverted. Then, as we 
discussed just above, the universe undergoes the following phases if we neglect particle 
production; 


— Phase 0 : Inhation takes place for 

TTleff < M 0 Mp (^^^n-2 J ^ (2.20) 

Due to the non-minimal kinetic term, the slow-roll conditions are easier to be 
satished than usual. Note that this condition is also the same as Mp <C H^/M. 

— Phase 1 : 0 starts to oscillate while the kinetic term is dominated by the non- 
minimal one for 


M < M < meff (2.21) 

In this regime, the energy density violently oscillates with time. When cj) 
dominates the universe, this means that the Hubble parameter also violently 
oscillates with time. Notice that mefr is an increasing function with time in this 
regime. 

— Phase 2 : The kinetic term of (j) in Eq. fl2.4p is dominated by the standard one for 

H <^M (2.22) 

In this regime, the evolution of the universe is the same as the ordinary case of 
the minimal oscillating scalar with Einstein gravity. 


Note that the Hubble parameter can have a large oscillating part in this case. In fact, 
the behavior of H can be classified into the following three regimes: 

( mesH for M < H (Phase 1), 

H ~ rUeaH^/M^ for MVmefr <H <M (Phase 2), (2.23) 

\ for H < M^/rriea (Phase 2). 


The reason why we consider these two cases will become clear in the next section. In the rest 
of the paper, we focus on the oscillation regime (Phases 1-2 for Case B) unless otherwise 
stated. 
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3 Perturbation 


3.1 Scalar perturbation 

In this subsection we derive the quadratic action for the scalar perturbation when the inflaton 
0 dominates the universe. The calculation basically follows those done in Refs. [3 El HU, 
although they focused on the inflationary regime. Below we will see that the perturbed 
quantities exhibit qualitatively different behavior in the inflaton oscillating regime. 

We must hrst solve the mixing between the scalar part of the metric and the inflaton. 
We use the Arnowitt-Deser-Misner (ADM) formalism, in which the metric is taken as 

ds^ = —N‘^dt^ + {dx^ + (d'^dt) [dx^ + (d^dt) , (3.1) 

where N is the lapse function, is the shift vector and is the 3-dimensional spatial 
metric. We take the unitary gauge 0(t, x) = (j){t) in the calculation. Then, the action is 
given by 


S = 






(3.2) 


In the ADM formalism, the Ricci scalar and the (OO)-component of the Einstein tensor are 
decomposed as 

R = ^ - E^) - {Kn^) - (3.3) 

^ , (3.4) 

where Eij is related to the extrinsic curvature Kij as 

By = /VA'y = 1 (vf'ft + vf 7 . - 7 y) , (3.5) 


the superscript “(3)” denotes covariant quantities with respect to the spatial metric 'jij and 
■jij is used to take the summation with respect to the space indices i,j ,.... Moreover, is 
the unit normal vector of the time-like hyper-surface, whose components are given by 


n‘ 




N 


(3,6) 


Thus, we rewrite the action as 


S = I dx 




N + 


2NMlM‘^ 



N 2N^M^M^ j ^ NMj 


(3.7) 
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Note that the action does not include time derivatives of the lapse function and the shift 
vector. Therefore the derivatives of the action with respect to them yield constraint equa¬ 
tions, which can be used to relate N and to the other scalar perturbations. The constraint 
equation for N is 









30^ ^ 

j 


0^ 


2V 

Jdl 


0, (3.8) 


and that for 0* is 



2N^MlM'^ 


(E'i - H'lE) 


= 0 . 


Now, we expand the metric 



(3.9) 


N — 1 “h Oil 

I3i = diip, 

-fij = a{tYe^‘^Sij, (3.10) 


where we have hxed the remaining one scalar degree of freedom by the gauge transformation. 
To zeroth order in perturbations, Eq. fl3.8p is nothing but the Friedmann equation Eq. 02.41). 
To hrst order, we obtain the following relations among a, 0 and ( from Eqs. 03.8p and Od.Op J^q 


where 


a = ^C. 


E = 


, F. , 02 F^G ■ 

^ = + 8,x = a—^C., 

iwi;' 


(3.11) 


(3.12) 


See Eq. 02.Op for the dehnition of e. By substituting these relations to the original action 
and performing some integration by parts, we obtain the following quadratic action for C, 
and hence the following discussion is not affected by this subtlety. 


#5 


S = MlM^ 


dx'^ a' 


,E^G 

~1p 


^ {d.cr 


(3.13) 


Since the vector and tensor perturbations do not couple to the scalar perturbations in the lowest order, 
we neglect them here. 

Although we write down the solution for ip here, ip does not contribute to the quadratic action because 
it appears linearly in the quadratic action after integration by parts. 

At the end-points of the oscillation, this action may be ill-defined since 0 = 0 at those points. However, 
the expression of the sound speed squared is the same in other gauges as we have shown in App. [Cj 




















where c* is the sound speed, which is given by 


c, = 


K = 


1 




+ 


3i72 




3F 


+ 


QH 

Jp 




^ 3H^ 1 + |e 

1 -I- — 

M2 1 - U 



(3.14) 

(3.15) 


where we have used Eq. fl2.7p to express 0 in terms of H. Note that due to Eq. fl2.6p . F‘^G/H‘^ 
in Eq. fl3.13p is always positive and no ghost mode exists in the scalar perturbation. For the 
sound speed squared c^, however, the situation is different. The term proportional to H/M'^ 
in the numerator in Eq. 03.141) can be positive or negative due to the violent oscillating 
feature of the Hubble parameter, and hence can be negative if this term dominates over 
the other terms. We consider Cases A and B separately. 

• Case A : In this case, one always has H < M and H ~ Then we estimate the 
sound speed squared as 


c 


2 

s 


1 + 0 



(3.16) 


where we have used 0^/MpM^ ~ for H < M. Therefore, the sound speed 

squared is positive definite and no violently growing mode exists in the scalar per¬ 
turbation. 


• Case B : From Eq. 02.23p . the sound speed squared in Phase 1 is given as 

~ ^ for M < H. (3.17) 

Even in Phase 2, the term proportional to H/M"^ dominates over the other terms for 
M{M/< H < M, and hence the sound speed squared is given as 


for M(M/meff)^'^^ < H < M. 


(3.18) 


Note that they are oscillating functions between positive and negative values in these 
regimes. This means that the sound speed squared becomes negative for some period 


during the oscillation]^ Hence, the scalar perturbation grows violently and the theory 


has an instability. In particular, the shortest wavelength modes are more likely to be 
enhanced, and hence we may need some UV completion in order to discuss what really 
happens in this case. Such a discussion is beyond the scope of this paper. Thus we do 
not argue whether this is really problematic or not in this paper. 


It is known that in some classes of the generalized Galileon theory the sound speed squared can be 
negative in the oscillation regime |15) . 
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In summary, an instability due to the negative sound speed squared exists in the scalar 
perturbation when the non-minimal derivative coupling dominates over the minimal kinetic 
term. If we require such an instability not to exist, the situation reduces to Case A where 
there is no Phase 1. This is why we consider both Cases A and B in this paper. 


3.2 Tensor perturbation 

Next, we consider the action for the graviton. We take the metric as 

= —dt'^ + a^{t) (e^) dx^dx^, (3.19) 

where hij is the transverse and traceless {dihij = ha = 0) part of the metric and is the 
matrix exponential of h. The relevant terms in the action are 




grav 


= J d^x^f^ 


Ml 


-R + 


2M2 


G 


00 


Since the Ricci scalar and the Einstein tensor are expanded as 


R = 6{2H^ + H) + - 


4 


hij j (X (^dihij^ 




) T ^ i^^ihij^ 


up to second order in gravitons, we get the following action: 

1 6^ 


'S'grav = f d'^X 


1 - 




hi,] -^{l + \TiT7^\(d,h„) 




2 M2 M2 


(3.20) 


(3.21) 

(3.22) 


(3.23) 


where we keep only quadratic terms in gravitons. Note that due to the non-minimal deriva¬ 
tive coupling, the graviton is non-minimally coupled to the oscillating scalar held. In par¬ 
ticular, the sound speed of the graviton is superluminal. However, due to Eq. fl2.6p . the 
graviton does not have the ghost-instability and the gradient-instability. 


4 Particle production 

In the present models, in particular in Case B, the Hubble parameter has a large oscillating 
part. In Ref. na, it is shown that the oscillation of the Hubble parameter causes production 
of scalar particles and the graviton. This is because all particles are coupled to the scale factor 
if they are not Weyl-invariant. Thus, in this section, we investigate particle production due 
to the oscillation of the Hubble parameter. Moreover, the non-minimal derivative coupling 
introduces a direct coupling between the inhaton 0 and the graviton. We also consider 
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the graviton production due to this direct coupling. We do not consider the production of 
fermions and vector bosons here because their couplings to the Hubble parameter vanish in 
the massless limit and hence are suppressed by small mass parameters. 

After summarizing general features on particle production, we discuss production of scalar 
particles and the graviton for both Cases A and B. We estimate the particle production rate 
assuming that the coherently oscillating inflaton (j){t) dominates the universe. However, 
as already mentioned. Case B has an instability associated with the negative sound speed 
squared which may invalidate this assumption. This fact should be kept in mind when 
considering the particle production for Case B. 


4.1 General discussion on particle production 


Before going to the discussion of the particle production in the present model, let us sum¬ 
marize basic properties of the particle production due to a coherently oscillating scalar field 
0 with the mass rries [TBl - lTB] . We call the produced (bosonic) particle as y here. In general, 
when a held couples to the coherently oscillating scalar held, the mass of the held also os¬ 
cillates because of that interaction. Due to the periodic dispersion relation, it is well known 
that the coupled held can be copiously produced [H], whose behavior can be seen in the 
stability/instability chart, for instance, of Mathieu equation for the sine or cosine oscillation. 
However, in the following, we will not dig into technical details of the equation under the 
periodic background. Instead, we demonstrate that its typical behavior can be understood 
more intuitively. 

We dehne a resonance parameter q as | Am^| = qmlg, where Am^ is an oscillating part of 
y’s mass squared due to the 0-dependence and |Am^| is its amplitude. Then, the interaction 
term can be written as^^ 




int 


Amlxl 




4>cxl, 


(4.1) 


where 0c and Xc are the canonically normalized helds, which may be diherent from the 


original 0 and y by an overall factor]^ In addition, <I)c is the amplitude of 0c. From this 


interaction, we can easily estimate the perturbative decay rate of 0 to y as 




87r<h2 


(4.2) 


Note that we have assumed g ^ 1: for g < 1, the “decay” of 0 into y particles is kinematically 
allowed, and hence the above estimation is justified. Otherwise, more elaborated treatment 

For interactions of the type with even powers of 0, such as A0^y^, it may be better to view the 
process as an annihilation, i.e. dint ~ /^c)4>cX^■ Th® estimation below, however, remains intact after 

interpreting the energy of y produced by the decay (annihilation) as meff/2 (mes). 

In the case with the non-minimal derivative coupling, 0c ^ H(j)/M for H/M ^ 1, and 4>c = 4> for 
H/M < 1. 
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of the particle production is needed [IS]. In this paper, we will only encounter the situation 
with g < 1. 


Up to here, our arguments do not depend on whether or not the resonance can take 
place. In the following, we focus on the narrow resonance regime (g ^ 1), which originates 
from the induced emission of previously produced y-particles. If the resonance occurs, it is 
well-known that there are resonance bands depending on the momentum k oi x [IB] • In the 
resonance bands, the wave function of y experiences exponential growth oc due to the 
coherent oscillation of (j). The growth rate /i can also be related to g as follows. The time 


evolution of the number density of y for t > l/(gmeff) is given as^^ 




(4.4) 


where and = p^/rries are the number density of y and (j), respectively. For the 
coherently oscillating scalar held 0, the number density is estimated as 




V 

rries 


meffd) 


(4.5) 


where we assumed V ~ If we dehne the momentum width of the resonance band 


as Ak, the distribution function can be estimated as 

2Tr‘^n^ 
Jx 


k‘^Ak 


(4.6) 


Since the band width arises from the oscillation of the mass Ak is estimated from 

kAk ~ g?^eff- Also, the momentum at the resonance band is given as fc ~ meff. Note that 
the momentum width from the mass oscillation must dominate over that from the decay 
width of (j), 


gmeff>r^, (4.7) 

for the resonance to take place Collecting all these ingredients and Eq. fl4.2p together, 
we can estimate the growth rate p from Eq. fl4.4p as 

fl 

/i = — ~ qmes. (4.8) 

At first, for -^1, the distribution function has relatively broad width of Akjk ^ l/(mefft) due to 
the uncertainty principle and grows as oc F, but the number density of y follows the ordinary evolution 
equation m, 

~ n^T^. (4.3) 

After t ^ 1/(grrieff), the band width becomes Ak ~ qnies. At the same time we have I, and after that, 
the induced emission becomes significant. Hence, we can use (14.41) for t > l/{qmes)- 
In our case, this assumption is justified both for H/M ^ 1 and H/M <C 1. 

We can prove this as follows. Assuming that (p completely decays into y and > qnies, the maximum 
phase space density of y is given by ^ n^/{k'^Ak) ~ ^ Thus we obtain 

j^max) ^ which means that the resonance does not happen for F^ > qm^s. 
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This implies that the number density of y exponentially grows as oc This is why we 
call q the resonance parameter from the hrst line. We would like to stress here again that 
technical details are not required to understand these properties. 

In the expanding universe, in order for the resonance to occur, the wave function must 
be enhanced enough before their momenta redshift out from the resonance band. Typ¬ 
ical timescale for the momentum within the band Ak escapes from the band is Atn ~ 
Ak/{kH) ~ q/H. Requiring that At^ is longer than the timescale of the resonant enhance¬ 
ment ~ l//i, we obtain [TB] 


> H. (4.9) 

It is noticeable that this condition ensures qm^g > ior q < 1 and < H. Eqs. fl4.2p . 
fl4.8p and fl4.9p are extensively used in the discussion in this section. 

4.1.1 Examples 

Here we briefly see the validity of our formula presented above. We take the potential to be 
V{4>) = Let us first consider the following coupling: C = with canonically 

normalized scalar held 0 and y. For this coupling, we easily obtain q = From 

Eq. 04.21) the decay rate of 0 is estimated to be ~ which is a well known 

result. 

Next let us consider the following coupling, 

C = -F\<l)){dxf. (4.10) 

where F{(j)) is an arbitrary function of <f>. In this case, by dehning Xc = F{(j))x, we obtain 
the following equation of motion of Xc 

fF p 

d\c-^Xc = Q. (4.11) 

In terms of Xci h seems that as if it obtains a mass term in the time-dependent background 
(/)(f). Thus, for homogeneous 0 held, we obtain q ~ F/{Fm‘^) ~ Hence, for the 

simple case of F = 1 0/A with |0| <C A, we obtain T^ ~ m^/(87rA^) from Eq. fl4.2p . 

reproducing the known results of perturbative decay rates. 

4.2 Production of scalar particle 

As we have seen in Sec. [21 in particular in Case B, the Hubble parameter has a large oscillating 
part in the present model. In Ref. [12], it is shown that 0 can decay into other non-Weyl 
invariant particles through the oscillation of the Hubble parameter since 0 is related to the 
Hubble parameter via the Friedmann equation. Thus, in this subsection we estimate the 
decay rate of 0 to another scalar particle y through the oscillation of the Hubble parameter. 
We take the mass of y to be negligible here. 
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The action of x field is given by 


= J ■ (4.12) 

We take only the background part of the metric as 

^x = \ j - a~‘^{dixf] . (4.13) 

Then, as done around Eq. 04.111) . we can deduce the resonance parameter as 

gW ^ JL^ (4.14) 

where we have used the fact that H is always larger than or comparable to H^. Thus we 
estimate the decay rate as 




rUeS^l ' 


(4.15) 


Now let us consider Cases A and B separately. Note that rries < H is always satisfied 

and hence parametric resonance does not occur for y. 


4.2.1 Case A 

First let us consider Case A. In this case, we always have H < M and the canonically 
normalized amplitude is the same as the original one: = $, and H ~ mes_^/Mp. In this 

case, we have 




Mj, 


(4.16) 


This is the same as production rate in the pure Einstein gravity limit, as it should be. To 
conclude, for Case A, the scalar production rate after inflation is always the same as that of 
the pure Einstein gravity case with minimal kinetic term. Here, note that even in the pure 
Einstein gravity, the Hubble parameter has a small oscillating part which causes the particle 
production [T^ . 


4.2.2 Case B 

Next let us consider Case B. In this case, at Phase 1 {H > M), the canonically normalized 
amplitude is given as $c = H^/M, and /MpM"^ ~ 0(1). At Phase 2 we have $c = $, 
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and H ~ nies^/Mp. Thus, by using Eq. fl2.23p . we obtain 

for M <H, 

for M’^/rries < H < M, (4-17) 

for H < M^/rrteff. 

Therefore it reduces to the pure Einstein gravity case with the minimal kinetic term for 
H < /rries- Note again that we estimate the production rate assuming that the coherent 
oscillation of (j) dominates the universe, although the instability of the scalar perturbation 
due to the negative sound speed squared may invalidate this assumption. 


T 




$2 

M4M|, 

Mj, 


4.3 Production of graviton 

Now let us focus on the graviton. Graviton production has two contributions: one is from 
the oscillation of the Hubble parameter, which results in the approximately same rate as 
those studied in the previous section, after reinterpreting x the graviton. The other 
contribution comes from the non-minimal derivative coupling term, in which the graviton 
and the inflaton directly couple with each other. Hereafter in this subsection we study the 
latter contribution. 

As we saw in the previous section, the action for the graviton is given by Eq. fl3.23p . 
From now on, we neglect the scale factor a in the action. This is because the amount of the 
graviton production due to the scale factor is expected to be the same as that of the scalar 
particle production, which we have already studied in the previous subsection. Thus, we 
consider the following action: 


^ — 
*-^grav 


1 

2 



1 G f 

2M2M2 j 


iG 


2M|,M2 


{dih 


■ijj 


(4.18) 


where we have rescaled the graviton as Mphij/2 —)■ hij. We further dehne another time t' as 


where 


dt = 


—dt', 

9 



1 02 
2 M2 M2 ’ 


2 _ 10^ 

^ ^ ^ ^ 2 M2 M2 ■ 


(4.19) 


(4.20) 

(4.21) 
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Note that p is positive definite (see Eq. fl2.6p h By using this time coordinate, the action is 
rewritten as 

' “ 2 


^ — 
*-^grav 


dt'dx^yjl 

2 V 4f‘ M},M‘ 


dhi- 


{dih 


ijj 


(4.22) 


Then, as done around Eq. 04.111) . we can deduce the resonance parameter g, 


(h) 


non-min 


as 




non-min ’ 


(4.23) 


Q: 


{h) 


where we have taken only the leading term of 6“^ /This treatment may be justified 
even in Phase 1 because of the relation 02.6p . From this gnon.^in given by 

1 for iP > M, 

(4.24) 

for H < M. 

Therefore ^ 1 is always satisfied. Note that, as already mentioned, the graviton 

production also occurs through the oscillation of the Hubble parameter in the same way as 
the minimally-coupled scalar field y. This contribution to the resonance parameter g for the 
graviton is denoted by g^;^ = g^^i. Thus an effective resonance parameter is given by 


(h) I (h) (h) 

g^ ^ = max { Cn-min 


(4.25) 


Since we have already seen g^^i < 1 in all cases, we also have g^^i ^ 1 in all cases. Then the 
decay rate of the infiaton into the graviton is given by 


T 


4>^h 


(g('l))' 


m 


eff 


$2 


(4.26) 


In the following, we consider the graviton production in Case A and B separately. 


4.3.1 Case A 

In this case, it is easily shown that 


Sh) 




{h) 

Q 

^inon-min 




> 1, 


(4.27) 


because H < M and mefr < M are always satisfied in this case. Therefore, in Case A, the 
graviton production rate in the infiaton oscillatory regime is always the same that in the 
pure Einstein gravity limit, as it should be 


T 


(p — 


Mj, ■ 


(4.28) 


Note also that () mes < 77 is also satisfied, hence no parametric resonance happens. 


In this sense, the difference between t and t' can be neglected. 
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4.3.2 Case B 


In Case B some non-trivial phenomena are likely to happen. In this case we have 


„{h) 


1 

El 

Ml 


m 


eff 


for M <H, 

for <H <M, 

for H < M‘^/mes- 


(4.29) 


In the following, we consider the graviton prodnction in Phase 1 and 2 separately. 

First, let us consider Phase 1. In this case, the non-minimal kinetic term dominates over 
the usual one, and hence the canonically normalized scalar held is related to the original 
one as <hc ~ H^/M. Moreover, note that ~ C>(1) is always kept in Phase 1 

when (j) dominates the universe; recall the discussion given around Eq. fl2.17p . Therefore, 
the perturbative decay rate is given as 




$2 • 


(4.30) 


Note that this is an increasing function with time. More importantly, Eq. fl2.17p means that 


7'“’ ~ 0(1). 


(4.31) 


is always kept in Phase 1 unless the amplitude $ is strongly damped due to the particle 
production. Therefore, the resonance condition fl4.9p is satished in Phase 1, and hence the 
resonance is expected to be induced as long as 0(t) can be regarded as homogeneous back¬ 
ground. In such a case, the perturbative decay picture is no longer useful, and a careful 
treatment is needed. It would lead to the exponential growth of the graviton with wavenum¬ 
ber around ~ rries. The production could be so efficient that universe may be dominated by 
gravitons, leading to the decrease of the inflaton amplitude <I> so that it saturates the reso¬ 
nance condition fl4.9p . These arguments are based on the assumption that (j){t) is regarded as 
homogeneous at least for the time scale of the growth of the graviton . How¬ 

ever, this assumption might be violated due to the instability associated with the negative 
sound speed squared. We do not go into details of this case. 

Next, let us consider Phase 2 with < H < M. In this case, the canonically 

normalized held (pc is identical to the original held (p. Therefore, the perturbative decay rate 
is given by 






(4.32) 


Note that is a rapidly decreasing function with time in Phase 2. Therefore, the resonance 
does not last long after entering Phase 2. 

Here is a comment. We obtained (x and T(j,^h oc (except for the $ dependence 
of nies) although the coupling between (p and the graviton in the original action fl3.23p is 
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of the form of from which one may naively expect oc and oc $2. This 

property may be understood in the language of the original action fl3.23p as follows. First, 
we must note that the coupling in Eq. 03.231) is in the following form: 


C 


int 


M2M2 


hij ) + {dihijY 


(4.33) 


If we consider decay of (j) particle at rest into two gravitons, the amplitude would be propor¬ 
tional to — |p|2 E and p are the energy and momentum of the produced gravitons, 

and we have used the 4-momentum conservation. For the usual on-shell gravitons, E"^ — |p|2 
would vanish since the graviton is massless. In the present case, however, the dispersion 
relation of the graviton is also modihed due to the non-minimal derivative coupling, and 
hence we get 


1 ^^ 


m 


eff 


$2 




(4.34) 


Therefore, the decay width is expected to be proportional to Again note that we simply 
neglect the effects of the instability here. 

Finally, for H < M’^/rries, the graviton production rate becomes the same as the pure 
Einstein limit fl4.28p . 


5 Conclusion 

In this paper, we have investigated the particle production in the inflationary models with the 
non-minimal derivative coupling between the inflaton and the Einstein tensor, 

focusing especially on the inflaton oscillation regime. In this model, when the non-minimal 
part dominates the inflaton kinetic term (Phase 1), the Hubble parameter violently oscillates 
with time if the inflaton oscillation dominates the universe. The oscillation of the Hubble 
parameter causes particle production, and we have estimated the decay rate of the inflaton 
in that process. More importantly, the non-minimal derivative coupling introduces a direct 
coupling between the inflaton and the graviton. We have seen that, via the non-minimal 
derivative coupling, the resonant graviton production may occur due to the inflaton oscilla¬ 
tion in Phase 1. As a result, the energy density of the inflaton may be efficiently transferred 
into that of the graviton, and the graviton may dominate the universe. 

However, in Phase 1, the sound speed squared of the inflaton rapidly oscillates between 
positive and negative values. It may result in the strong instability for the shortest wave¬ 
length mode. It may or may not be phenomenologically problematic, but analysis in such a 
violent situation is beyond the scope of this paper. A conclusion is that if we simply avoid 
such an instability within this effective held theory, the particle production rates are the 
same as those in the models with Einstein gravity with minimal kinetic term as found in 
Ref. [12]. Also it should be noticed that inhation models with non-minimal derivative cou¬ 
pling taking dominant role such as the new Higgs inhation model [7] necessarily experience 
this instability after inhation. 
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Although we have assumed in this paper that the non-minimal coupling between the 
inflaton and the Einstein tensor takes the form C ~ the graviton production is 

expected to occur in more general G^^-type models (Gs-terms in the context of the gener¬ 
alized G-inflation 0 )- This is because a direct coupling between the inflaton and graviton 
is introduced in these models, and hence the graviton is expected to be produced efficiently, 
as in the case of the present model. However, we should also be careful on the instability of 
perturbations in these models. 
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A Inflation with non-minimal derivative coupling 


In this Appendix we briefly discuss slow-roll inflation when the non-minimal derivative part 
of the kinetic terms takes a dominant role. Sin ce G °° ~ and G®-^ ~ —3H^/a? during 


inflation, the Lagrangian of inflaton 0 is given aq^lf 

377^ 


C ~ 


2M2 


{d4>) —V{(j)) ~ 


A0" 




(90)2-R(0). 


We can define canonically normalized field 




y/nMMf 


(n+2)/2 


and the Lagrangian is rewritten as 


(A.l) 


(A.2) 


/:^--(90)2-R(0), R(0)= - 


A\ 


n 


MMf 


2n 

n+2 


(A.3) 


Note that this transformation is meaningful only in the slow-roll regime. After inflation, we 
cannot approximate as 3MpH^ ~ V, hence such a transformation is not possible. This is 
the difference from the running kinetic inflation [8]. 

Precisely speaking, ^ and hence the sound speed of S(j) is modified. But its effect is small 

as long as we are concerning physical quantities up to the first order in slow-roll parameters. 
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Therefore, in terms of 0, the potential just behaves like oc The slow-roll 

equation of motion is 

+ P = 0, q ^ (A.4) 

According to the standard procedure, slow-roll parameters at the e-folding number N is 
calculated as 


1 

2' 




n 

V) 2(n + 2)iV’ 


W = = 


n — 2 


V 2(n + 2)N 


(A.5) 


The slow-roll conditions evi\'r]v\ ^ 1 are equivalent to rries <C H. We obtain the scalar 
spectral index and the tensor-to-scalar ratio as 


Us = 1 — 6ev + 2r]v = 1 


2(n + l) 
{n + 2)N' 


r = 16ev = 


8n 


[n + 2)N' 


(A.6) 


These are inside the 2a range of the Planck constraint for n < 4 [3l[Tll|20]. In order to ht 
the CMB normalization = 2.2 x 10“® observed by Planck [3], we need 


m^M 

Ml 


~ 2 X 10"^° 



for n 


2 


(A.7) 


and 

~ 1 X 10-3^ for n = 4. (A.8) 

Note that in order for inflation to take place in the non-minimal regime, we need the 
following condition: 


3P2 


am"-2 

M2 


2 

n+2 


> 1 , 


(A.9) 


where 0e is the held value at the end of inhation. This implies S> M for n = 2 and 
A S> {M/MpY for n = 4. Thus we can take S> 10^^ GeV for n = 2 and A S> 10“^° for 
n = 4 by choosing M appropriately while keeping the CMB normalization consistent with 
observations. 


B Proof of adiabaticity 


In this Appendix, we prove the adiabaticity of J. 


The energy conservation Eq. fl2.10p reads 


d 

dt 




+ 3M 



3M d 
JpJt 



0 . 


(B.l) 
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The second term is of the order of ~ 0{Hp^), and hence we safely neglect this term here. 
Thus, we must hnd a new adiabatic invariant from the equation 


A 

dt 




3H d 
JpJt 



~ 0 . 


(B.2) 


In order to hnd an adiabatic invariant in this system, we dehne the following quantity /; 


I = 


1 




(B.3) 


where Ci and C 2 are some numbers we determine below. By taking derivative with respect 
to time, we get 


dt dt 


NP 2 


9 — C 2 H d 
2 




(ciC2 — 9ci — C 2 + 9) — 3MpCiHH, 


(B.4) 


where we have used 0^/2 + V = 3MpH^{l — 3(fr/2MpM^). Therefore, for 

Cl = 1, 

C2 = 3, 

the following equation is obtained: 


(B.5) 

(B.6) 


— (/ + 3MlH) ~ 0. 


(B.7) 


This means that the quantity I + 3MpH is adiabatically conserved. This quantity is related 
to J as 


2J = I + 3MlH. (B.8) 

Therefore, we have proved that J is an adiabatic invariant. Note that the adiabaticity of J 
holds for all Phases 0-2. A more systematic way to derive an adiabatic invariant for general 
actions is described in m 

C Sound speed squcired in another gauge 

In this appendix, we show that the formula of the sound speed squared of the scalar pertur¬ 
bation fl3.14p does not change in another gauge which is well dehned even at the end points 
of the oscillation, 0 = 0. We again use the ADM formalism, where the metric is given as 

ds^ = —N'^dt^ + '^ij {dx^ 13^dt) {dx^ + jd^dt) . (C.l) 
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Here we take the gauge condition as 


N — 1 -f- 
A = diii, 

'yij (i(t') 

(p = 0(t) + 5(p, (C.2) 


where we have neglected the vector and tensor modes since they do not affect the quadratic 
action of the scalar perturbation. Here (p{t) is the background part of the scalar field in this 
gauge. Note that we can take this gauge condition even at 0 = 0 points. We treat a^ip and 
6(p as perturbations and expand the action with respect to them. Then, the action is given 
to the second order as follows: 


S= (Tx 


1 3H^'. 

2 (1 + igr I'*■#' 


9H'^ 

-<P{l + j^]a6cP- 


2a2 
2H(p 


1 + 


2H 

+ 


[dMf - 




M2 M2 
diadi6(j) — V'a6(j) 


a 


a2 M2 


^^(50 - 2HM^p ( 1 - - f , \a\d’^'ip 
M2 ^ 2 M2 m2 ' ' ^ 


(C.3) 


where we have omitted the background part of the action. Also, first order terms vanish due 
to the background equation of motion. 

Two important features can be read from the above action. Firstly, a and ip does not 
have the kinetic term, and hence their equations of motion give constraint equations among 
6(p, a and ip. Secondly, the action depends only linearly on ip. This means that, after solving 
the equation of motion, ip does not contribute to the original action. Thus, we need not 
explicitly solve the equation of motion for a to get the form of ip in terms of 6(p. This feature 
greatly simplifies the calculation. 

By solving the equation of motion for ip, we obtain the following relation; 


0 

2MM2 (1 - |e) 



H - 


(C.4) 


where e is defined in Eq. 02.61) . After substituting this into Eq. PC.31) and doing some 
integration by parts, we obtain the following kinetic term for 5(p\ 


Fkin = m2 m2 


f, 4 

/ dx a ^—— 

<502 - 4 {dMf 

/ 02 

1 Ci \ 


(C.5) 
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where the definitions of F, G and are the same as those given in Sec. 13.11 and we have 
omitted the mass term for 5(j) since it is not relevant to the sonnd speed sqnared. In order 
to obtain this expression, the following background equation 


F / _ 1 A 3F2\ H . 

M^V 27 2 V ^ M2 J ^ 2M2^’ 


(C.6) 


is useful. Note that the action flC.Sp is well defined even at 0 = 0 points since G oc Thus, 
we conclude that the formula of the sound speed squared is the same in the other gauge 
condition flC.2p . Note that, after including the mass term, the actions 03.131) and OC.Sp are 
equivalent once we identify ( and <50 as 

C = (C.7) 

0 


This is nothing but the gauge transformation from Eq. 0C.2p to Eq. 03.10p . 
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